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Abstract 


The purpose of this program was to study direct optical 
modulation techniques applicable to integrated optical data 
preprocessors. After reviewing several techniques potentially 
useful in direct modulators, it was decided to concentrate our 
research effort on the analysis and fabrication of a field effect 
type modulator. A series of computer modeling studies were 
performed to determine the effects of semiconductor claddings on 
the fields of propagating waves in planar dielectric waveguides. 
These studies predicted that changes in the propagation characte- 
risitics of waveguides clad with silicon and gallium arsenide 
could be made large enough to be useful in modulators. These 
effects are dependent on the complex permittivity and thickness 
of the cladding. Since the conductivity of the cladding can be 
changed by the photon generation of hole-electron pairs, 
incoherent light may be used as the input to modulate a coherent 
light beam. Waveguides were fabricated and silicon claddings 
were applied to verify the theoretical predictions. 
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I . INTRODUCTION 


An integrated optics holographic comparator has been 
proposed (1) for use as a data preprocessor for airborne optical 
sensors. This device, as presently proposed, requires the 
conversion of the sensor optical signal into an electrical signal 
which then modulates a coherent optical beam by means of a set of 
metal electrodes deposited on a lithium niobate waveguide (Figure 
1). This technique introduces noise and extra hardv/are because 
of the needed signal conversion from optical to electrical. 

The objective of this investigation was to determine the 
feasibility of direct modulation of a coherent optical beam by an 
incoherent optical source to eliminate some of the above 
problems. Several direct modulation techniques were investigated 
during the first year of this study. One technique investigated 
was the metal-barrier-metal, or more commonly metal-oxide-metal 
device (MOM), which makes use of tunneling effects through a 
barrier region. While the MOM device shows great promise for 
this application, it is still in the early stages of research and 
its ability as an optical modulator in an integrated optical 
circuit has yet to be investigated. Two other techniques 
currently being developed are the Pockels-Readout-Optical 
Modulator (PROM) type memory and display devices and a spatial 
light modulator being developed by M.I.T. These devices were 
discussed in Section 2 of the first annual report (2) . 

The concept selected for investigation during the latter 
part of this study was that of a photon induced change in conduc- 
tivity in a semiconductor clad optical waveguide (Figure 2). The 
background investigation and preliminary computer prediction on 
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Figure 1. Schematic of One Configuration of the Integrated Optics 
Holographic Comparator. (Ref. 1) 
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Figure 2. Photoconductive Modulator 
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the feasibility of using a clad waveguide type modulator are 
presented in detail in references 2 and 3. Section 2 of this 
report briefly summarizes the results of these previous studies 
along with attempts at experimental verification of the predicted 
photon modulation effect. 

Section 3 of this report discusses the results of computer 
predictions on the effects of changing the real part of the 
complex permittivity of the semiconductor cladding. It was 
orginally postulated that deposited silicon films being used in 
the experiments might have a significantly different permittivity 
than that used for the computer predictions. The predicted 
results did not verify this postulate and thus additional 
experimental studies have been initiated as discussed in Section 
7. 

In order to better understand the results of the computer 
predictions and to develop a theory for the predicted modulation 
effects, the fields in the waveguide, substrate and cladding 
region were calculated and plotted. Section 4 presents the 
results of these calculations and the multilayer interference 
theory used to explain the computer predictions. 

Unfortunately the modulation scheme being studied suffered 
from an extremely high attenuation due to the semiconductor 
cladding behaving as lossy dielectric. In Section 5 the results 
of attempts to reduce this attenuation through the addition of a 
buffer layer between the waveguide and the cladding are 
presented. The computer predictions indicate a significant 
reduction is possible through the use of a proper buffer layer. 
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Because of a number of problems which were experienced in 
depositing good photoconductive silicon films, selemium-clad 
waveguides were investigated. It is shown in Section 6 that 
selemium willproduce the same photon modulation effects as 
predicted for silicon and can thus be used if a suitable deposi- 
tion systems is found. 

Section 7 discusses experimental results of numerous studies 
on the temperature degradation occuring in ion exchanged wave- 
guides and the design of new experimental measurement devices. 
The effects of exposing ion-exchanged waveguides to elevated 
temperatures for short periods of time were investigated to 
determine what methods of depositing semiconductor films on the 
waveguides were feasible or whether sputtered waveguides would be 
required. The design of a new series of experimental devices to 
be used in future measurements is also discussed in this section. 
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II. SEMICONDUCTOR-CLAD DIELECTRIC WAVEGUIDES: SUMMARY OF 
PREVIOUS RESEARCH 

In this research, a planar waveguide configuration has beejn 
used for all computer modeling and fabrication. A semi-infinite 
glass substrate (refractive index Nj) covered first with a wave- 
guide material of refractive index N2 such that N2 > then a 
layer of semiconducting material forms the waveguide structure 
(Figure 2 ). The order of the dielectric and semiconductor 
layers is reversible. The fourth layer is air in all cases and 
for all calculations the waveguide material was polystyrene al- 
though any dielectric with a similar refractive index will 
produce the same results. 

This basic structure has been analyzed utilizing classical 
electromagnetic theory and has been found to support propagation 
in the waveguiding region with a set of fields charcterized by; 



where y = ct + j 3, propagation constant 

u) = radian frequency of propagation wave 
a = attenuation constant (nepers/meter) 

e = phase constant (radian/meter) 

Each propagating path is characterized by its own 
propagation constants and In order to understand a 
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structure's physical behavior, we must solve for a and B. The 
relationship we must solve is the dispersion relationship, which 
relates a and 3 to the physical properties of the structure and 
to the characteristics of the exciting wave (i.e. frequency). 
For the previously described structure, such a relationship has 
been derived (4). The general technique and result for this 
structure is as follows (2) : 

(1) It is assumed that all materials except the 
semiconductor are approximately lossless and that 
propagation is in the z direction. 

(2) The nonzero components of the TE and TM modes are 
determined and substituted into the wave equation. 

(3) The wave function is obtained by solving the wave 
equation. By invoking the condition of continuity 
of the wave function and by forcing the field 
components to satisfy Maxwell’s boundary condi- 
tions, the dispersion relation is obtained; 


"-X3 *=3 = 



j t -1 *^34 ^x4 , 

+ tan — + niT 


K 


x3 


( 2 . 2 ) 


where K^j = i for TE modes 

Kij = Ej^/ej for TM modes 
Pxi = - €[ K^) + j(eV 2a3)]^^^ 

Kxi = [(e. - 3^ + a^) + j (2a3 - eV K^) 

n = 0, 1, 2 = mode order 
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Kq = 2tx/\ 

£i =s ej^' - j£j[" = complex permittivity of region i. 

Because this is a transcendental equation, it has been solvable 
only by numerical techniques. Two separate approaches have been 
used to solve this, and both have given identical results. For 
the remainder of this work, primary attention will be focused on 
a and g for the guiding region (N3) • 

One attempt to solve the relation involves a razor search 
technique used by Handler and Macdonald ( 5 ), In this program, 
the operator selects initial guesses for a, 8 and convergence 
factors for error estimation. The computer then uses a random 
walk technique to calculate values for a and B satisfying the 
dispersion relation. Results are based on input of the above 
mentioned data plus the physical structure parameters set at the 
beginning of the run. This program (called WAVE) has one inhe- 
rent weakness due to the requirement for finding the tangent 
function. While the tangent may be found by at least two diffe- 
rent methods both produce questionable results as points of 
discontinuity are approached or when the tangent function changes 
quadrants. This was a problem of concern when t2 became very 
small (t2 << 1 yni) . 

A second method of solving for a and B is to formulate the 
problem as an eigenvalue problem as described by Smith and 
Mitchell (6) . In this approach, characteristic matrices are 
created for each layer based on given material parameters. These 
matrices are then solved for the required eigenvalues and thus 
a and B. The solutions obtained with this program (called 
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MODEIG) are identical to those obtained with WAVE at all points 
plus it is able to compute points at the tangent discontinuities 
which WAVE cannot. 

Before beginning v/ork on this grant, changes in the attenua- 
tion of progagating modes in a semiconductor clad waveguide 
(Figure 2) had been observed as the semiconductor cladding 
material was varied (Figure 3). Any lossy material has a complex 
permittivity given by 

g = £'..j£" = £»_-] 4rra 

0 ) 

where e' is the real part of the permittivity and a is the 
conductivity of the material at frequency co. This can also be 
expressed in terms of the complex refractive index 


where k is the extinction coefficient. 

It was noted that the real part of the complex permittivty 
is about the same for germanium (17.4) and gallium arsenide 
(14.3), but the conductivity at A = 632.8 nm is quite different 
[Ge, o = 3.2 X 1q 2 (ohm cm)“^ and GaAs, a = 25.4 (ohm cm)”"^]. 
Figure 3 shows the predicted variation in attenuation for these 

two materials when the cladding thickness is semi-infinite (t 2 “ “ 
in Figure 2b). While this was not a largechange in a, it 
indicated an important phenomenon to be investigated. The fact 
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Figure 3. Attenuation of Ge and GaAs Clad Waveguides 
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that a change in the conductivity of the cladding material 
changesthe attenuation of the propagating made in the waveguide 
then leads one to predict modulation effects since the 
conductivity of a semiconductor can be changed through photon 
generation of hole-electron pairs. This has recently been 
demonstrated in microwave guides using semiconductor cladding 
(7). 

Having noted the effect of a change in conductivity on the 
attenuation of a waveguide, the next step was to try to enhance 
the magnitude of the effect through proper choice of cladding 
materials, cladding thickness and waveguide thickness (t 3 ). 
Allowing t 3 to vary as the conductivity was changed only shifted 
the characteristic attenuation curves (as in Figure 3) up or down 
and produced no detectable change in the phase characterisitc of 
the waveguide. 

Previous experience indicated that the attenuation normally 
decreased as the cladding thickness (t 2 ^ was decreased so the 
next step was to vary the cladding thickness while also varying 
the conductivity to see if this had any beneficial results. The 
results obtained were perplexing since they were entirely 
different than anything previously observed. Figure 4 shows the 
results for the GaAs clad waveguide as t 2 is varied. Not only 
does the attenuation vary with conductivity but it also seems to 
oscillate after the cladding thickness decreases past 1.0 ym. 
This oscillation had never been observed before so it raised a 
number of questions about possible errors in the computer 
program. The most significant result of these predictions v/as 
that the conductivity ppeared to have a relatively large affect 
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Figure 4. Variation of Attenuation with GaAs Cladding Thickness 
for the TEq Mode (t 3 = 1.0 um) 


12 



on the attenuation at certain cladding thicknesses. Since only a 
few thicknesses were used in calculating Figure 4, additional 
points were calculated to detemine the oscillation frequency. 

Figures 5 and 6 show the variation of attenuation with t 2 
for a fixed conductivity (normal value) after the additional 
points were calculated. Now the oscillation is even more 
pronounced and there seems to be a discontinuity in the mode 
index curve at 0.09 ym. At thicknesses greater than 0.1 ym the 
mode index curve appears to be continuous. As noted before, this 
oscillatory behavior had never been observed before in any of our 
calculations but after verification using both computer programs 
the conductivity (a) and cladding material were varied. The 
three most common semiconductors were selected to see what effect 
the material permittivity had on the results. Since GaAs had 
already been investigated, silicon and germanium were used as 
claddings. The complex permittivities and refractive indices are 
shown in Table I which follows. 

Table I 

Semiconductor Parameters at A = 632.8 nm 


Material Permittivity Refractive iBdgx, 



£' 

e" 

n 

k 

GaAs 

14.3 

1.21 

3.79 

0.16 

Silicon* 

16.76 

1.75 

4.1 

0.213 

Germanium* 

14.43 

19.54 

4.4 

2.22 


*values for amorphous thin films 
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Figure 5, GaAs/Air-Clad Waveguide Attenuation vs. t^ 

(TEq mode, normal conductivity; t^ = 0.8 pm) 
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Note that the relative magnitudes of the real and imagninary 
parts of the complex permittivity of gallium arsenide and silicon 
are the same. Germanium, however, has a considerably higher 
imaginary part (conductivity) than either silicon or gallium 
arsenide. If material conductivity has any effect on the charac- 
teristics observed it should be evident in germanium. Figures 7, 
8, 9, and 10 show the results for silicon and germanium. As 
might be expected, silicon exhibits characteristics nearly iden- 
tical to those of GaAs. Germanium shows almost none of the 
oscillations which were characterisi tcs of GaAs which can be 
attributed to the large conductivity used for germanium in this 
calculation. 

The next step in the investigation was to vary the conducti- 
vity of one of the materials and plot the change in attenuation 
and phase as a function of cladding thickness. Figures 11 and 12 
are the results for silicon where the percentage change in 
attenuation is plotted with the percentage change in conductivity 
as a parameter. The curve indicates that for a proper choice of 
silicon thickness ( ,13 microns) the attenuation changes by 
approximately the same percentage as the conductivity. The phase 
characteristic (Figure 12) also shows an appreciable change for 
proper choice of silicon thickness. It is also interesting to 
note that maximum phase change occurs at the 0% crossing of the 
attenuation curve, which indicates a device could be designed as 
either an amplitude or phase modulator. 

Based on these calculations several devices were constructed 
to verify the predicted attenuation variation with cladding 
thickness. It was hoped that these devices would not only verify 
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Figure 7. Silicon/Air Clad Guide Attenuation vs. t^ 
(TE^ Mode, normal conductivity) 
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the calculations but prove that the magnitude of the effect was 
large enough for use as a modulator. Silicon was selected as the 
semiconductor to be deposited on the waveguide because of the 
well -developed deposition technology. Discussions with several 
people working with amorphous semiconductors indicated potential 
variation in the optical properties of these amorphous films with 
deposition method and temperature. Thus values of refractive 
index and loss coefficients used in our calculations may not be 
accurate for real films. After initial investigation of the two 
principal methods of deposition of amorphous silicon (CVD and e- 
beam evaporation) it was decided to use the e-beam technique 
because of better reproducibility in optical constant. The e- 
beam system at Cornell University, Materials Science Department, 
was used for our initial films. 

The film thickness was selected based on the predicted 
attenuation and phase variations. A wedge shaped film varying in 
thickness from 0.1 to 1.0 micrometers was deposited with the 
Cornell system in the configuration shown in Figure 13. A 
detailed discussion of the deposition technique is available in 
reference 3. According to the predictions, the attenuation 
should vary by >900 dB over this thickness region and thus the 
effect would be readily observable. 

This device was fabricated by first constructing an ion 
exchange waveguide (3) and then evaporating the silicon film on 
the waveguide through a mask. The silicon thickness gradient was 
obtained by evaporating through a v/edge shaped oscillating 
shutter. De'ktak measurements of the resulting films indi cated 
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Figure 13. Experimental Device for Verification 
of Oscillations in Intensity. 
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a film thickness varying from .06 to .105 micrometers rather than 
the 0«1 to 1.0 desired. Measurement of the attenuation of the 
best of three guides fabricated was accomplished by coupling a 
laser beam into the waveguide through a prism coupler. The 
intensity of the guided beam was then measured with a photo- 
detector and a liquid coupler before and after the silicon film 
giving the relative attenuation of the silicon film. 

Due to the silicon film not being in the thickness range 
desired and problems with making repeatable measurement using the 
liquid coupling technique, the measured data v/as difficult to 
correlate with the predicted curves. Figure 14 shows the measured 
points obtained from two seperate sets of measurements of the 
guide. There appears to be a general decrease in the attenuation 
as the thickness approaches 0,105 micrometers but it is not a 
well defined curve. The predicted data indicated the attenuation 
should be increasing in this region. It was felt that there were 
two possible causes for the discrepancies between the measure- 
ments and the predictions. The first explanation v/as that an 

oxide layer of 300A had grown on the silicon thus reducing the 

o o 

film to (a thickness of) 300A (to 750A). This would make the 

experimental curve match the calculated curve in magnitude and 

shape. Calculation of oxide growth on bulk silicon indicates a 
o 

10-20A layer should be present. Unless the thin film behaves 

differently from the bulk material or if there were an interac- 

o 

tion at the glass silicon interface this 300A change in 
thickness needed to match the data with predictions could not be 
justified. The second explanation was that the silicon film was 
contaminated by impurities in the silicon and in the vacuum 


25 



Attenuation 

(dB) 



4 - Run 1 
X Run 2 



% Average 

Rst . trend 



Figure 14. Experimental oc Measurements 




chamber, thus changing the real part of the permittivity. A 
theoretical investigation of this second postulated effect is 
discussed in detail in section 3. 

The work reported in the remainder of this report occurred 
from October 1980 to June 1981 and is an attempt to explain the 
origin of the effect. This work includes theor eticalanalysi s, 
additional computer predictions and attempts at fabrication of 
improved experimental devices. 
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III. EFFECT OF CHANGES IN CLADDING PERMITTIVITY 


Several explanations for the discrepency between experimen- 
tal results and computer predictions were proposed. One possibi- 
lity was that some changes had occurred in the real part of the 
permittivity of the silicon film from the value used in the 
computer calculations. Possible explanations for the 
permittivity change include mechanical irregularities in the Si 
film and inadvertent doping of the semiconductor layer. 

Since the mechanisms for the e* shifts were not well under- 
stood, a trial-and-er ror method was used to determine the value 
necessary to correlate experimental results with computer calcu- 
lations. Initially, a Ae' = 16.76 - 15.76 = 1.0 was used; 
however, the magnitude of the curve shift was insignificant. 
Further attempts were made and the results of the computer calcu- 
lations are shown in Figure 15 . Observe that the real part of 
the permittivity must change by at least 40% in order to produce 
an appreciable horizonal curve shift; hov/ever, the attenuation- 
thickness curve shifts vertically as well (the local minima 
shifts but the local maxima remains the same). In addition, a 

40% decrease in the permittivity (from 16.76 to 10.76) only 

o o 

shifts the curve 200A, unfortunately, a 500A curve shift is 
necessary to achieve any kind of correlation between experimental 
results and computer calculations. This indicates that a permit- 
tivity decrease on the order of 80% would be necessary to achieve 
correlations; but, a permittivity decrease also shifts the magni- 
tude of the local minima, and further correlation problems are 
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introduced. It appears, then, that an e* shift is not a suitable 
explanation for the discrepency between experimental results and 
computer predictions. 
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IV. ANALYSIS OF PREDICTED WAVEGUIDE CHARACTERISTICS 
A. Theory 

The damped oscillations of the attenuation-silicon thickness 
curves and the phase-silicon thickness curves can be explained in 
terms of both a ray-optics approach for light propagation in the 
v/aveguide and a careful analysis of the electric and magnetic 
fields of the propagating wave within the structure. 

Heavens (8) considered a thin film of material with n = 2 
and k (k =aA/4Tr) varying and calculated the expected phase change 
on reflection of an incident wave at an air film surface. For k 
very small, the phase change oscillates, and for k large, the 
phase changes very little until the film thickness approaches 
zero. This behavior has also been predicted by Stratton (9) and 
is similar to the effect we have observed. 

The proposed explanation for the results is that there isan 
interference effect across the waveguide and when constructive 
interference occurs in the semiconductor cladding, the 
propagating wave field strength will be greater and higher 
attenuation is expected. Likewise, destructive interference will 
cause a much lower attenuation. Normalized phase shifts must be 
at a maximum when the interference effects are at a minimum, 
similar to the results of transmission calculations of waves 
through an absorbing medium. Increasing the conductivity of the 
semiconductor cladding' increases the field strength attenuation 
in the cladding which decreases the effect of the interference. 
At the points where a greater conductivity causes a lower 
attenuation, the greater conductivity increases the reflection at 
the cladding-core interface, and thus there is less field 
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strength in the absorbing medium and less attenuation in the 
propagating v;ave. 

Accurate electric and magnetic field plots of the wave 
within the structure were obtained by using program MODEIG (6) 
and confirm the validity of this argument. Program MODEIG ob- 
tains a 2 X 2 characteristic matrix of Maxvi/ell's equations and 
boundary conditions for each waveguide layer defined as: 


cos + (i Y^)sin 
Yjj)sin oos 


where 




(the complex phase 
thickness of the layer) 



(the complex propagation constant 
of the layer) 


Y ^ = generalized parameter of the field 
distribution within the layer 


Through a complex series of multipications, of the above charac- 
teristic matrix for each layer, the propagation constant andthe 
incremental field distributions within the structure are calcu- 
lated. From the characteristic matrix, we observe that the field 
distribution in each waveguide layer is a complicated 
trigonometric function of the complex propagation constant. 

For a. lossless structure, however, the wave function profile for 
each layer is effectively reduced from its complicated matrix 


32 



representation to the following: 

j 

e"^4 i = 4 i 

A- cos (a^ + y) i = 3 ' 

, _ j (ojt-K z) 

^ A2Cos(a2 + 6)i = 2 | 

Aj e“l 1 = 1 ; 

We again note a trigonometric variation in the semiconductor 

cladding and dielectric guide, and a decaying exponential in the 

semi-infinite layers. 

B. Calculated Field Distributions 

We now consider the problem of correlating the local 

maximum/minimum points on the attenuation-silicon thickness curve 

with the electric, and magnetic field distributions at the 

appropriate thickness. Results from program MODEIG indicate 

o 

that the presence of a thin silicon film (<100A) has little 
effect on the wave function profiles; the profiles are very 
similar to those of the three-layer structure. For the thick 
silicon film structure, however, the lowest order mode of the 
lossless three-layer structure appears to couple to the modes 
assosciated with the semiconductor-dielectric interface. 
Furthermore, the coupling between the modes supported by the 
three-layer lossless structure and the modes associated with the 
semiconductor-dielectric interface determines the attenuation of 
the complete four-layer structure. 


33 



The real part of the wave function profile in the transverse 
direction for the TEq mode is shown for a cladding thickness, 
t 2 =.007 ym (Figure 16). We note that the thin silicon film does 
not appreciably affect the wave function profile of the TEq mode. 

For the first local minimum (t 2 = .05ym), we observe that 
the field strength at the silicon-dielectrc interface approaches 
zero (Figure 17). We also note the expontially decaying solu- 
tions in the outer, semi-infinite layers. 

For t 2 = .09 ym (the first local maximum), we observe a sharp 
peak in the wave function profile at the silicon-dielectric 
interface (Figure 18). 

Similar behavior is noted for the next local minimum/maximum 
pair (t 2 = .13 ym and t 2 = .18ym) (Figure 19 and Figure 20). The 
field strength is effectively zero at the silicon-dielectric 
interface for the local minimum, and a sharp peak in the v/ave 
function at the interface is evident for the local maximum. We 
also note an increase in the number of field oscillations in the 
semiconductor cladding. 

Again, for the next local minimum/maximum pair (t 2 = .22ym 
and t 2 .26ym) (Figure 21 and Figure 22) we observe an increase 
in the number of field oscillations in the silicon and a corres- 
ponding increase or reduction in field strength at the silicon- 
dielectric interface. Interestingly enough, at the local maximum 
(t 2 = .26pm) the field in the dielectric guide decays rapidly to 
zero, and the field strength is at a maximum at the silicon- 
dielectric interface. This behavior v/as not evident for the 
other local maxima (t 2 = .09ym and t 2 = .18ym). A considerable 
field strength was, of course, apparent at the interface; how- 
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ever, a sizable field was still noted in the dielectric guide. 
Consequently, fields were calculated at two thicknesses (t 2 “ 
0.255ym and t 2 = 265vim) on either side of t 2 = 0.26iim. For 
these two thicknesses we observe that there is, again, a sizable 
field contribution in the dielectric guide (Figure 23 and Figure 
24). This result, in particular, indicates that the local 
"maxima" used for the computer calculations are not, precisely 
local maxima. Time and monetary constraints prevented more 
accurate determination of the local maxima points; however, 
appropriate field distributions for these thicknesses are still 
apparent. 

The field plots indicate that the attenuation of the four- 
layer structure may be explained quite simply as a coupling 
between the basic three-layered waveguide (air-guide-substrate) 
modes and the surface wave mode supported by the silicon- 
dielectric interface. For a thin silicon cladding (.007ym) the 
wave function is not appreciably distorted and is everywhere 
positive which is appropriate behavior for the TE^ mode. For a 
local minimum on the attenuation-silicon thickness curve, the 
field at the semiconductor-dielectric interface is zero, and for 
a local maximum, a considerable field is set-up at the interface 
which decays rapidly through the guiding region and the semi- 
infinite substrate. Furthermore, we note that the number of 
field oscillations in the silicon cladding is related to the 
silicon thickness and consequently to the order of the local 
maxima/minima points (first, second, third, etc.). 
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We recall that the wave function profiles in layers 2 and 3 
are trigonometric for the lossless case. Indeed, the structure 
of interest is not lossless, so a slight distortion from the 

idealized trigonometric profiles appears. By solving for the TEq 
mode of the four-layer structure we are requiring one half-cycle 
for the wave function profile in the guide; however, by 
progressively increasing the silicon thickness we are obtaining 
more oscillations for the trigonometric profile variation within 
the semiconductor. If the silicon thickness is such that the 
profile can go through exactly one cycle (or an integer multiple) 
we have perfect coupling at the silicon-dielectric interface and 
a minimum of the attenuation-silicon thickness curve; however, if 
the thickness is such that the wave function profile goes through 
half-cycle multiples, we observe a maximum on the attenuation- 
thickness curve. We note, of course, a slight distortion from 
the pure trigonometric variation in the silicon because boundary 
conditions must be matched at the appropriate interfaces. 
Furthermore, we observe that the direction of the field (positive 
or negative) in the dielectric is dependent on the tangent of the 
profile in the semiconductor (see t 2 = .255, .260, .265 ym). 

Finally, we note that the period of the attenuation-thickness 
curve (.08ym) is approximately identical to the half-period of 
the sinusoidal wave variation in the silicon cladding. A slight 
discrepenecy occurs, of course, as the silicon thickness is 
varied and the propagation constant, t <2 =w2 ye^ _ (K^/k^)^ for the 
silicon layer changes slightly; consequently, the period of the 
oscillations of the attenuation curve is not always 8 nanometers. 
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The effect that we have described is quite similar to cer- 
tain developments in the theory of grating couplers (6, 10, 11, 
12, 13). These results indicate that for some configurations 
(the grating layer placed between a surface guide and a high 
index substrate), the coupling efficiency is very dependent upon 
the grating period and guide thickness. For some combinations of 
parameters, the coupling efficiency could be large, but for other 
slightly changed values, the efficiency could become essentially 
zero. In reference (6) the authors state that these effects can 
be described as periodic coupling betv/een the guided mode and 
other leaky modes of the same guide. Thus, if Bji is the propaga- 
tion constant for the bound mode, and 3 2 real part of the 
propagation constant for a leaky mode then strong coupling to the 
radiated beam (from the leaky mode) can be expected when the 
phase matching condition is satisfied, K = 2 tt/A = 132^-321, where 
A is the grating period. 

In summary then, our results indicate that the damped oscil- 
lations of the attenuation-silicon thickness curve can be ex- 
plained both in terms of a ray-optics approach for light propaga- 
tion and in terms of a careful analysis of the electric and 
magnetic fields within the waveguide. Both approaches yield 
identical conclusions. 
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V. REDUCTION IN ATTENUATION: FIVE-LAYER STRUCTURES 


Thin dielectric buffer layers have been used to lov/er the 
attenuation losses of metal-clad dielectric waveguides (14). 
These layers are placed between the dielectric core and the 
metal, and act as buffers to remove a large portion of the field 
from the metal claddings. 

The effect of an S 1 O 2 buffer layer on the attenuation - Si 
thickness characteristics will be considered. Where an Si02 layer 
has been added between the guiding region and the semiconductor 
cladding to reduce the attenuation. 

Using program MODEIG, five different cases were studied. 
They are; 

Case i. verification of known results (convergence to 

three and four layer cases for ^Si 02 = °° and 
^Si02 

ii. (a) vs Si thickness for different Si 02 
thicknesses (e r ” 2.12) 

iii. (a) vs Si thickness for different Si02 
thicknesses (£r = 2.387) 

iv. (a) vs Si thickness for different Si02 

0 

pemittivities (Si02 thickness = lOOOA) 

V. (a) vs Si thickness for different Si02 

0 

permittivities (Si02 thickness = 2000A) 
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Case i 


The attenuation constant (a) vs Si thickness is plotted for 
an Si02 (^r “ 2.12) thickness of lOOA (Figure 25). The 
attenuation-si thickness curve for the simple four-layer struc- 
ture is shown for reference. Note that the attenuation-Si 
thickness curve for an Si02 thickness of lOOA converges to the 
simple four-layer case, which is, at least, intuitively correct. 
Further computer mode], ing studies indicate that the 

attenuation-si thickness curve converges to the lossless three- 

0 

layer case for ^si02 ^ 2000A and that the attenuation constant 

(a) can be effectively reduced to zero with large buffer layer 
thicknesses. 

Case ii 

The attenuation constant (a) vs Si thickness is plotted for 

0 

two different Si02 (^r = 2.12) thicknesses (1000 and 2000A) 
Figure 26, The attenuation-Si thickness curve for the simple 
four-layer structure is shov^^n for reference. The familiar damped 
sinusoidal behavior is apparent for all buffer layer thicknesses 
(100, 500, 1000 and 2000A). Most importantly, the attenuation can 
be reduced significantly with the addition of the Si 02 buffer 
layer while still preserving the damped sinusoidal behavior of 
the attenuation-Si thickness curve. The percentage decrease in 
attenuation (relative to the four-layer values) for different 
buffer layer thicknesses is summarized below. 
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Figure 25. Attenuation Characteristics for 5-layer Structure 
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Figure 26. Attenuation Characteristics for 5-layer Structure 


•50 



cont . case ii 


Si Thickness (um). ^2 = 2.12) thickness. 




0 

lOOA 

0 

500A 

0 

lOOOA 

0 

2000A 

.05 

(local min) 

-2.5% 

-17% 

-40% 

-7 5% 

.09 

max 

“5.8% 

-30% 

-56% 

-84% 

.13 

min 

-2.4% 

-17% 

-3 9% 

-7 5% 

.18 

max 

-1.5% 

-14% 

-37% 

-74% 

.22 

min 

-2.8% 

-18% 

-41% 

-7 6% 

.26 

max 

-2.6% 

-18% 

-42% 

-76% 


A significant decrease in attenuation is obtained by increasing 
the Si 02 thickness. 

Case iii 

The attenuation constant (a) vs Si thickness is plotted for 

0 

two different Si02 <£r = 2.387) thicknesses (1000, 2000A) (Figure 
27). The attenuation-si thickness curve for the simple four- 
layer: structure is shown for reference. The familiar damped 

sinusoidal behavior is apparent for both buffer layer thicknesses 

'0 

(1000, 2000A), and the attenuation can be reduced significantly 
by a buffer layer. The percentage decrease in attenuation, 
relative to the four-layer values, for different buffer layer 
thicknesses is summarized below. 
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SILICON THICKNESS tg (/xm) 


Figure 27. 


Attenuation Characteristics for 5-layer Structure 







cont . case iii 


Si 

Thickness (\im) 

S iOo (^ r) — 2 . 387). 

thickness 



0 

0 



lOOOA 

2000A 

05 

(local min) 

-25% 

-48% 

09 

max 

-30% 

-54% 

13 

min 

-25% 

-48% 

18 

max 

-23% 

-47% 

,22 

min 

-2 5% 

-49% 

,26 

max 

-25% 

-49% 


Again, a significant improvement is obtained by increasing the 
Si 02 thickness. 

Case iv 

The attenuation constant (a) vs Si thickness is plotted for 
different Si02 (thickness = lOOOA) permittivities (£ r ~ 2 . 12 ., 
2,387) (Figure 28). The attenuation-Si thickness curve for the 
simple four-layer case is shov;n for reference. The familiar 
damped sinusiodal behavior is apparent for both cases, and the 
attenuation can be reduced significantly. 

The percentage decrease in attenuation (relative to the 
four-layer values) for the different permittivities is summarized 
below. 
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ATTENUATION a (n/m) 



Figure 28. Attenuation Characteristics for 5-layer Structure 


o 

iii..IhiffKnesg £1^2 -(thickness = lOOOA) Permittivity. 

= 2.12 ~ 2.387 


05 

(local min) 

-40% 

-25% 

09 

max 

-56% 

-30% 

13 

min 

-39% 

-25% 

18 

max 

-37% 

-23% 

,22 

min 

-41% 

-25% 

,26 

max 

-42% 

-25% 


Si 02 (£r “ 2.12) gives a more significant attenuation decrease. 

Case V 

The attenuation constant (“) vs Si thickness is plotted for 

0 

different Si02 (thickness = 2000A) permittivities (e r = 2.12, 
2.387) (Figure 29). The attenuation vs Si thickness curve for 
the simple four-layer case is shown for reference. The familiar 
damped sinusoidal behavior is apparent for both buffer layers, 
and the attenuation can be reduced significantly. 

The percentage decrease in attenuation (relative to the 
four-layer values) for the different permittivities is summarized 
below. 
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SILICON THICKNESS tg (/xm) 


Figure 29. Attenuation Characteristics for 5-layer Structure 




2.387 



.05 

(local min) 

-7 5% 

-48% 

.09 

max 

-84% 

-54% 

.13 

min 

-7 5% 

-48% 

.18 

max 

-74% 

-47% 

.22 

min 

-76% 

-49% 

.26 

max 

-76% 

-49% 


Again, Si02 (er 2.12) gives a more significant attenuation 
decrease. 

These computer modeling studies indicate that the attenua- 
tion can be reduced significantly with the addition of an Si 02 
buffer layer between the guiding region and the semiconductor 
cladding. The familiar damped sinusoidal behavior is still 
evident for all five-layer cases considered; however, a slight 
decrease in the relative amplitude of the oscillations is 
apparent as the buffer layer thickness is increased which indi- 
cates a slight decrease in the attenuation shift for a given 
conductivity change. We also note that the attenuation can be 
more effectively reduced with a lov/er permittivity buffer layer. 
This is, at least, intuitively correct since the index of refrac- 
tion of the buffer layer is less than that of the guide, and a 
lower index implies more confinement of the radiation to the 
guide. Conversely, computer studies indicate that the attenua- 
tion increases if the buffer layer has a higher index of 
refraction than the guide. Studies also indicate that a buffer 
layer increases 3/K slightly, but decreases the amplitude of the 
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oscillations on the 6/K-silicon thickness curve. For a practical 

device, then, an overall attenuation reduction with Si02 buffer 
layer must be carefully weighed against a decrease in the 
attenuation shift for a given conductivity change. 
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VI. SELENIUM-CLAD WAVEGUIDES 


Several other semiconductors exhibit a damped sinusoidal 
behavior as the cladding layer thickness is increased. One 
substance, elemental selenium, appeared ideal for device applica- 
tions. 

Amorphous silicon, in contrast to the chalcogenide glasses 
(i.e. selenium), contains a large number of dangling bonds which 
form intergap states and fast recombination centers. The 
recombination sites limit the conductivity changes for all but 
the brightest intensities. One solution is to "passivate” with 
hydrogen which bonds to the silicon and effectively reduces the 
number of recombination centers (15) . 

Diffraction studies of selenium, hov/ever, indicate that the 
glass is comprised largely of helical chains and the evaporated 
or sputtered thin films have varying proportions of rings, as 
opposed to dangling bands; hence, passivation is not necessary. 
Furthermore, thin film deposition of selenium (melting tempera- 
ture = 200*^0 can be easily accomplished by quenching from the 
melt (16, 17). 

The preparation independent optical properties of selenium 
have been determined. At A = .6328 ym, the complex permittivity 
t = 7.9 + jl.O. Computer modeling of a selenium-clad dielectric 
waveguide again revealed the familiar damped sinusidal behavior 
of the attenuation vs thickness curve (Figure 30). Selenium, then, 
appeared ideal for device applications. 

Unfortunately, further investigation revealed that thin film 
(<.3 ym) preparation from the melt is impossible and that this 
technique creates toxic selenium gas. 
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Figure 30. Attenuation Characteristics for Selenium-Clad Waveguide 





VII. EXPERIMENTAL STUDIES 

A. Ion-Exchange Waveguide Temperature Degradation 
Several possible semiconductor deposition techniques were 
investigated including RF sputtering and vacuum deposition. 
Unfortunately, many methods involve subjecting the ion-exchange 
waveguides to high temperatures. Of particular interest, then, 
was the temperature degradation of the ion-exchange guides. 

Fabrication procedures have been described in (18) and v;ill 
only be briefly reviev/ed. The basic technique for ion-exchange 
optical waveguide fabrication is based upon the diffusion of 
silver ions into soda-lime glass and the replacement of sodium 
ions which diffuse out of the glass. The index of refraction 
(traditionally, a microscope slide) is raised, thereby forming a 
waveguide. The number of propagating modes is proportional to 
the diffusion time. 

Several guides were prepared, and a table relating the 
number of propagating modes with the diffusion time is shown 
below. 
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Waveauide 

Number of Modes 

Diffusion Time 

1 a and 

t> 

4/4 

4 

min 

2 a and 

b 

4/4 

4 

min 

3 a and 

b 

5/5 

6 

min 

4 a and 

b 

5/5 

6 

min 

5 a and 

t» 

1/1 

15 

sec 

6 a and 

b 

1/1 

15 

sec 

7 a and 

b 

1/1 

45 

sec 

8 a and 

b 

1/1 

45 

sec 

9 a and 

b 

7/7 

10 

min 

LO a and 

fc) 

7/7 

10 

min 


Some ol: these guides were heated to 300°C and then to 500®C, 
and the number of modes was observed. 


Label 

Modes (27°C) 


Modes (5Q0°C). 

Wave 1 a/b 

4/4 

4/4 

D 




E 

3 a/b 

5/5 

5/5 

S 




T 

5 a/b 

1/1 

1/1 

R 




0 

7 a/b 

1/1 

1/1 

y 




E 

9 a/b 

7/7 

7/7 

D 


We note that all waveguides were able to support the same number 
of modes even after heating to 300°C; however, all guides were 
effectively "destroyed" after heating to 500°C, The conclusion, 
then, is that RF sputtering still is the ideal semiconductor 
deposition technique. 
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B, Design of Experimental Devices 

Previous attempt to confirm the damp sinusoidal behavior of 
the four-layer silicon-clad optical waveguides have not been 
successful. Fortunately, an RF sputtering systems along with a 
silicon target was available for semiconductor deposition. 

Radio frequency (RF) sputtering is a technique in which a 
source material (a target) is bombarded with gas ions accelerated 
by a high voltage. Particles from the target are ejected, 
traverse the vacuum chamber, and are subsequently deposited on a 
substrate as a thin film. 

It was decided that a large number of waveguides would be 
sputtered with uniform silicon films ranging in thickness 
from .04ym to .40 ym in .Olym steps in order to confirm the damped 
sinusoidal behavior of the attenuation-silicon thickness curve. 
The uniform thickness silicon bar allov/s more accurate attenua- 
tion measurements because several measurements may be performed 
at different places on the guide and the results averaged, A 5 
inch diameter silicon wafer was soldered to an aluminum holder 
and sputtering began. 

A sputtering holder was designed and constructed which would 
allow the deposition of a silicon film lym wide and 20ym long 
on the waveguide surface. 

In addition, a Hall probe mask was built in hopes that 
several sputtered silicon material characteristics (the material 
type, dopant concentration, and mobility) could be determined. 
Furthermore, Hall mask dimensions were designed so that a quick 
and accurate resistivity measurement could be made using a four- 
point probe (19). 
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Unfortunately, amorphous silicon (sputtered) is permeated 
with dangling bonds which render the sputtered silicon photoelec- 
trically and photoconduct ively dead. This is less than ideal 
characteristics for an optical modulator which relies on photon- 
induced conductivity changes in the semiconductor cladding. 
These dangling bonds lie betv;een the valence and conduction bands 
and act as fast, non-radiative recombination sites. One solu- 
tion, hov/ever, is to sputter in hydrogen which effectively passi- 
vates the dangling bonds thereby revealing the interesting 
properties of silicon. Other atoms, including oxygen, nitrogen, 
and flourine may aid the passivation process. 

The available system, however, only allows RF sputtering in 
argon and oxygen. Argon is completely ineffective for passiva- 
tion, and oxygen is less than ideal. In fact, a number of 
sources expressed serious doubts about the possibility of passi- 
vation with oxygen. Fortunately, motivation v/as found in a 
number of obscure journal articles. 

Paesler, Anderson et al (20) present the results of dark 
conductivity, photo-conductivity, photoluminescence, and optical 
absorption experiments on films of sputtered oxygenated amorphous 
silicon-hydrogen alloys. The authors conclude that oxygenated 
a “ - Hjj films exhibit photoconductive and photoluminescent 
response characteristic of glow-discharge films (i.e. passivation 
occurs) . 

Shimizu, Kumeda et al (21) prepared amorphous silicon with a 
low number of dangling bonds without the incorporation of 
hydrogen. The systems was pumped dov/n to about 6 x 10“^ torr. 
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and 99.995% pure O 2 was introduced into the chamber to a pressure 
of approximately 1 x 10“^ torr. Argon v/as then introduced to a 
total pressure of 2 x 10”2 to 3 x 10“^ Torr. 

Ishii, Naoe et al (22) significantly reduced the number of 
gap states in amouphous silicon by sputtering in an argon-oxygen 
mixture (O 2 partial pressure « 1 x 10“^ torr). Photoconductive 
properties comparable to those of hydrogenated amorphous silicon 
v;ere observed. 

The currently available vacuum system can realistically 
only pump down to a base pressure of 1 x 10“^ torr, which is not 
sufficient to duplicate the previous experiments; however, sili- 
con v;as sputtered in an oxygen atmosphere (p = 1 x 10“2 torr) and 
a uniform layer of Si02 was obtained. Growth rate v;as 
approximately 1200A/hr. A Dektac measurement along with a color 
calibration confirmed the film's uniformity. Many other films 

were sputtered in Ar -02 mixtures. Growth rates of approximately 

0 

1500A/hr were obtained; however, tests indicated that the 
sputtered silicon was still photoelectrically dead so passivation 
was not realized. Finally, experiments were halted when the Si 
v/afer continued to come unglued from the holder. 

A new silicon sputtering target has been obtained and 
sputtering will begin again. Only confirmation of the damped 
sinusoidal behavior will be initially attempted (purely a bulk 
effect). Passivation, then, is not necessary, and argon will be 
used as the gas. 
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VIII. CONCLUSIONS 

This study has led to several important conclusions about 
the development of direct modulation devices employing 
semiconductor clad optical waveguides. Based on this study and 
the following conclusions we believe such modulator are 
practical: 

(1) Calculations using two different methods of 
solving the clad waveguide problem produce the 
same oscillatory behavior. 

(2) The percent change in attenuation or phase shift 
is large. 

(3) Similar experiments by other workers indicate the 
effect is real and useable in a device. 

(4) Attenuation can be reduced significantly by 

inclusion of a buffer layer. 

As pointed out in section 2 all calculations using both 
available computer programs verify the oscillating behavior of 
the attenuation and phase constant of the guided wave when clad 
by a thin silicon or GaAs films. The fact that the oscillations 
in phase and amplitude are 90° out of phase allows the device to 
be designed as either a phase or amplitude modulator. The 
particular application originally invisioned for the direct 
optical modulator would require a phase modulator however the 
amplitude modulator might find application in communication 
circuits. 
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The percentage change in attenuation and phase with 
conductivity was large enough to make such a device feasible. 
The following table briefly summarizes these results for silicon 
films. 


Film 

Thickness 

(urn) 

% Chance 

in 

Changes in 
(%) 

Phase Shift 
( rad/mm) 

.14 


50 


42 

<0.1 

.14 


25 


22 

<0.1 

.17 


50 


<5 

-4.8 

.17 


25 


<3 

-3 

.22 


50 


35 

<0.2 

.22 


25 


19 

<0.1 

A main concern 

is whether 

a film 

<0.25 micromet 

,ers thick will 

absorb € 

: n 0 u g h 

incident 

light 

to change the 

conductivity 

signif icantly. 

Previous 

calculations(23) have 

shown that the 

minimum 

predic 

ted change 

in condi 

activity for a £ 

silicon film is 

6.7% with 

, a 50% 

change in 

light intensity. Even 

this produces a 


0,6 radian/mm phase change and has provided encouragement to 
continue with experimental devices. 

Lee and his co-workers (7) have constructed an optically 
con trolled millimeter-wave switch using conductivity changes in 
dielectric-clad waveguides and successfully demonstrated the 
modulation effect. Because of the lower frequency of the pro- 
pagating wave, the silicon film was considerably thicker than 
that required for a direct optical modulator and thus a stronger 
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effect was observed than might be expected in thin films. From 
this data it appears that an optical device based on the same 
principal is worth building. 

The buffer layer studies in section 5 indicte the potential 
reduction in attenuation possible through the use of a thin 
silicon dioxide film between the waveguide and the semiconductor. 

A 2000A SiC )2 film reduces the attenuation by 7 5% so that a 1 mm 
long film (.,13 mm thick, Si) would have an attenuation of about 
2.1 dB which would be an acceptable loss for a phase modulator. 
In general the buffer layer allows us to shift the attenuation 
curve down while having no significant effect on the phase curve. 
This is an ideal control variable for either a phase or amplitude 
modulator. 

Finally it should be mentioned that experimental data 
obtained from the one device fabricated v;as not conclusive and 
unfortunately technical problems with Si sputtering prevent 
fabrication of additional devices. Any future work must include 
experimental devices to verify the predicted oscillatory 
behavior. It is also evident that since the effect we are trying 
to make use of is predominate at film thickness of less than 1 
micrometer the problem of absorption of the modulating signal 
must be addressed in greater detail. It may be necessary to 
couple the modulating signal into the thin semiconductor film 
through the use of a second guiding layer on top of the silicon 
film. These problems will be investigated in future work in this 
area. 
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